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Biquadratic antisymmetric exchange terms of the form — [Cijefj(si x Sj) 
vector connecting sites i and j and a — x,y, due partially to magnetoelectric coupling effects, are 
shown to be responsible for the spin-flop helical phase in CuFeC>2 at low magnetic field and tem- 
perature. Usual biquadratic symmetric exchange, likely due to magnetoelastic coupling, is found to 
support the stability of axial magnetic states at higher fields in this nearly-Heisenberg like stacked 
triangular antiferromagnet. A model Hamiltonian which also includes substantial interplane and 
higher-neighbor intraplane exchange interactions, reproduces the unique series of observed com- 
mensurate and incommensurate periodicity phases with increasing applied magnetic field in this 
highly frustrated system. The magnetic field-temperature phase diagram is discussed in terms of a 
Landau-type free energy. 

PACS numbers: 75.30.kz, 75.50.Ee, 75.40.Cx, 62.20.Dc 
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I. INTRODUCTION 

The unusual magnetic properties of the stacked tri- 
angular antiferromagnet (STAF) CuFeC>2 have been the 
subject of considerable investigation and speculation over 
the past several decades, especially since the recent dis- 
covery of an unconventional magnetoelectric (ME) cou- 
pling induced by relatively weak magnetic field pi** The 
focus of the present work is on the surprising sequence 
of magnetic states which occur at low temperature with 
increasing field strength applied along the hexagonal c- 
axisfi^^ These can be characterized with the follow- 
ing spin polarization vectors S and in-plane periodici- 
ties: linearly polarized S||c period-4 (P4), helically po- 
larized Sic incommensurate (IC), linearly polarized 
S||c period-5 (P5), axial and canted S period-3 (P3) 
phases. Although the rhombohedral R3m crystal sym- 
metry present at temperatures above the onset of lin- 
early polarized S||c IC magnetic order in zero field at 
Tjvi=14 K allows for axial magnetic anisotropy, usual 
spin-orbit coupling is absent in the S=5/2, L=0 magnetic 
state of the Fe 3+ ionsi 7 -^ (The extent to which Hund's 
rules apply to this semiconductor with non-trivial elec- 
tronic structure^ may be questioned.) Anisotropy may 
thus be expected to be small so that the origin of stability 
of the axial phases is not obvious. The existence of such 
a wide variety of magnetic states with different periodici- 
ties in this compound is in contrast with the other weakly 
axial STAF's such as CsNiCbj which show only period-3 
phases and a typical spin-flop transition^ Neutron scat- 
tering experiments also indicate a period-two modula- 
tion of the magnetic order in CuFe02 along the c-axis of 
the corresponding hexagonal unit cell suggesting substan- 
tial antiferromagnetic (AF) interplane coupling*^ Such 
a complex phase diagram defies explanation through a 
model Hamiltonian based on conventional magnetic in- 
teractions. 

Remarkably, some important aspects of the magnetic 



ordering (e.g., stability of P4 and P5 states) have been 
accounted for by a model 2D triangular Ising AF with 
very large second and third-neighbor exchange interac- 
tions, J 2 /Ji=0.45 and J3/ Ji=0.75, respectively, where 
Ji is the nearest-neighbor (NN) exchange i 11 ^ 2 ' 13 The 
physical origin of these unusually long-ranged exchange 
effects is unknown; 7 -^ but there is evidence from per- 
turbed angular correlation measurements that J2 and 
J3 are substantial^ Analysis of the model proposed 
here demonstrates the importance of inter-plane cou- 
pling, J'j 7 -^ for stabilizing the P4 state with smaller val- 
ues of J2 and J3, an effect which has been found in other 
frustrated systemsJ^ 



It has also been emphasized that spin-lattice cou- 
pling is an important effect in CuFeC>2, as observed 
in magnetostriction^ and in the pressure dependence of 
T/vjii Such interactions are also believed to be relevant 
for the structural phase transition to monoclinic C2/m 
symmetry which occurs as the temperature is lowered to 
about 11 K. It is near this temperature Tjv2 that a dis- 
continuous magnetic phase transition from the linear IC 
phase to the ground state P4 ordered structure occurs^ 
It is of interest to note that although the type of mag- 
netic order is quite different, a similar C/2m - R3m struc- 
tural phase transition in solid O2 has been shown to be 
driven by magnetoelastic coupling.- 9 Magnetoelastic cou- 
pling is known to give rise to biquadratic symmetric ex- 
change terms of the form Hq = — J2ij G( r ij)[ s i ' s j] 2 j 
where = rj — rjjSS and has been related to the P4- 
IC phase transition in Ho. 21 This type of coupling can 
also arise from higher-order effects within the Hubbard 
model— Such coupling terms favor linearly polarized 
magnetic order— and are shown here to be partly re- 
sponsible for the stability of the modulated axial phases 
with s||H in CuFe02 at high fields even in the absence 
strong anisotropy. 
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II. ANTISYMMETRIC EXCHANGE 

Key to understanding the magnetic phase diagram of 
CuFeC-2 is the occurrence of an induced electric polar- 
ization P in the IC spin-flop phased This type of ME 
effect^! is atypical due to the fact that the delafossite 
crystal space group R3m contains inversion symmetry. 
Coupling between noncollinear spins at sites i and j and 
the electric polarization can be shown to be compatible 
with R3m symmetry if it also involves the lattice vector 
Vij and is of the for m 25 i 26 

n c = \j2 C ( r »)(Pij x ey) ' (Sj x Sj ) z (1) 

ij 

where = r^/r^, C(— r) = C(r) and the subscript z 
indicates the z component, parallel to the hexagonal c 
axis. Such coupling involves only components P _L z and 
should be present in all crystal systems which contain a 
symmetry plane where the spin Hamiltonian is isotropic 
in s at second order. It has also been used as a mechanism 
for electric- field induced chirality selection. 2 ' 

Following Katsura et alM- and Bergman et al. 2S , con- 
sider now adding to the Hamiltonian the lowest order 
contribution in P as a NN sum Hp = \Ap J2<ij> Pfj- 
Elimination of P by minimizing Hcp = Hc + Hp results 
in a new type of biquadratic antisymmetric exchange in- 
teraction of the form 

n CP = -g^-E E [cfriK-fa >< *i) (2) 

P a <ij> 

where a — x,y. It is important to note that, as in the 
case of biquadratic symmetric exchange, these antisym- 
metric terms of the form ~ [(sj x s,-) • z] 2 can be deduced 
purely from symmetry arguments applied to a Hamilto- 
nian with only spin degrees of freedom^ In the case of 
CuFe02, ME coupling provides for one likely microscopic 
mechanism. The impact of He p on the evolution of spin 
states with increasing applied field is demonstrated be- 
low. 



III. MODEL HAMILTONIAN 

The full model Hamiltonian used here to describe the 
sequence ordered states of CuFe02 incorporates all of the 
effects described above. Weak axial anisotropy is intro- 
duced in the form of the exchange (two-site) type J z (r y ) 
with possible physical origins from Fe 3+ defects^ or cer- 
tain forms of antisymmetric superexchange exchange. 29 
Anisotropy is also included in the biquadratic exchange 
term. For simplicity, anisotropy of the single-ion type is 
omitted from the present analysis. With these consider- 
ations the Hamiltonian can be written as 

H = H., + H C p +H g + H Jz + H Gz - H • ]T s(r,) (3) 



where the first three terms are given above, bilinear ex- 
change anisotropy is of the form Hj z = J2ij Jz( r ij) s i s j 
and biquadratic exchange anisotropy is given by Hg z = 
— Y^ij G z (rij)[sf Sj] 2 . The last term is the Zeeman cou- 
pling to a magnetic field H taken here to be applied along 
the hexagonal c axis. In addition to the four quadratic 
Heisenberg exchange interactions J\, J2, J3, and J', NN 
in-plane and inter-plane magnetoelectric coupling Gr- 
and C , as well as NN in-plane and inter-plane bi- 
quadratic exchange couplings G\ and G' are included in 
the model calculations. The impact of modifications to 
the in-plane exchange coupling due to the monoclinic lat- 
tice distortion at low temperature is accounted for in a 
model proposed in Ref. [5|. This effect was explored 
within the present model and was not found to substan- 
tially alter the qualitative results presented above. 

As a prelude to a discussion of the specific model 
Hamiltonian, and complimentary to the results of Mekata 
et al^, the effect of interplane exchange coupling on 
stabilizing the P4 ground state phase is demonstrated 
here. For this purpose, the Ising model is adequate. The 
quadratic exchange contributions to the Hamiltonian are 
written as Hj = Y^ij J( r ij) s i ' s j where J(?ij) > for 
AF coupling and the sum is over sites within, as well as 
between, the ABC stacked triangular layers which form 
the hexagonal crystal representation of the rhombohc- 
dral structure. Ground-state phases of Ising spins s||z||c 
with non-zero J\ > 0, J2,^3 > and J' > were de- 
termined analytically assuming the spin structures de- 
scribed by Mekata et al^ and verified numerically using 
the local-field method of Walker and Walstedt^S Figure 
1 demonstrates that long-period P4 and P8 spin configu- 
rations are stabilized by interplane exchange even in the 
absence of J3 but that third-neighbor in-plane coupling 
does serve to enhance the stability of these structures. 
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FIG. 1: Ground state phases of the Ising model on a rhombo- 
hedral AF Ji > illustrating the effect of interplane coupling 
J' > on stabilizing longer-period structures described by 
Mekata et al.. 11 Solid curves correspond to J3 = and bro- 
ken curves J3 > 0. 

With only weak axial anisotropy, an applied magnetic 
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field tends to destabilize a linear AF spin configuration 
with sf | |z. An illustration of the effect of biquadratic 
symmetric exchange to enhance the stability of axially 
ordered states is shown in Fig. 2. Using the model de- 
scribed above the evolution (S z ) 2 = (1/A) ^-(.Sf) 2 (an 
average over all N sites) in the P4 phase with increasing 
H shows clearly this enhancement. 
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FIG. 2: Effect of biquadratic exchange on the axial component 
of spins in the P4 state averaged over all sites with increasing 
magnetic field. Model Hamiltonian (3) was used with 3% 
anisotropy and Ji = 1, J% = 0.33, J3 = 0.3, J' = 0.4, Ci = 
C" = 0. 



IV. APPLICATION TO CuFeOa 

The set of parameters required to reproduce qualita- 
tively the sequence of field-induced spin configurations 
observed in CuFe02 is not unique. Some guidance can be 
found by estimation of exchange and anisotropy strengths 
based on the saturation and spin-flop fields, respectively, 
as done by Petrenko et ali^ Using their approach and 
accounting for the six NN in-plane interactions, leads 
to Ji ~ 3K, with anisotropy being only a few percent 
of this value. For simplicity, all of the anisotropy con- 
stants were set to be 3% of their corresponding isotropic 
coupling strengths, i.e., Jkz — 0.03Jfe(k=l,2,3), J' z — 
0mJ',G u = 0.03Gi,G' z = 0.03G". All parameters 
are normalized to the NN in-plane exchange by setting 
J\ = 1. Other parameters of the model were chosen us- 
ing guidance from previous model results and to yield a 
semi-quantitative representation of the unique sequence 
of field- induced phase transitions observed in CuFe02- 

Ground state phases which minimize the full energy 
expression (3) were calculated numerically using the 
method of Walker and Walstedt 30 with ABC stacked tri- 
angular layers of dimensions L X L X M , where M = 1 
corresponds to a unit cell containing all three layers. Pe- 
riodic boundary conditions were imposed and spin con- 
figurations with in-plane periodicities ranging from L=l - 
L=18 and up to M=4 unit cells along the c axis were con- 
sidered. Absolute energy minima were found by compar- 



ing results among calculations using up to 50,000 differ- 
ent random initial spin configurations. With this finite- 
size method, in-plane IC states characterized by (q, q) 
modulations can be not be modeled directly. However, 
such states can be deduced by demonstrating a slight re- 
duction in energy as L increases if q ~ n/L where n is 
an integer. For example, if the model energy is a mini- 
mum for an IC phase with q ~ 4.5, the method will give 
a lower energy with L = 9 (n = 2) than either L = 4 or 
L = 5. 

Sets of parameter values which give the correct series 
of phases were determined partially through trial and 
error. One such choice is J2 = 0.33, J3 = 0.3, J' = 
0.4, Ci = 0.3, C = 0.1, A P = l,Gi = 0.06, G' = 0.02. 
The ground state energies as a function of field corre- 
sponding to selected lower-energy periodicities (L, M) 
are shown in Fig. 3 for this set of values. The lowest en- 
ergy states as H increases follows the observed sequence 
(4,2)-(IC,2)-(5,2)-(3,l), with transitions at approximate 
values of H c i ~ 0.75, H c2 ~ 1.75 and H c3 ~ 3.0, following 
the notation of Ref.[j|. Saturation is found to occur at 
H c5 ~ 13.5. The phase labeled IC was deduced to be in- 
commensurate since in this field regime L=4, 5, 9 and 14 
states have nearly identical energy, with L=9 being the 
lowest. This conclusion is consistent with the observed 
value of q ~ 2/9^. Numerical accuracy of the method is 
limited at larger values of L. Note that the L=4, 9 and 
14 phases have a simple periodicity of 2 along the c-axis 
(q z = 1/2, representing six triangular layers). Additional 
Fourier components q z are also present in the [5,2] state. 
These characterizations are consistent with observed neu- 
tron diffraction d&t&££& At higher field values, the [3,1] 
phase is found with M = 1 (only three layers are re- 
quired), followed by saturated ferromagnetic $f||c order. 
Although substantial values of longer-range exchange in- 
teractions Ji and J3 are required to yield the correct se- 
quence of phases, they are significantly reduced by inter- 
plane coupling from those assumed in Refs. 14]fl2lfl3| . 

Figure 4 shows the site averaged (S z ) 2 with increasing 
H associated with the phases which minimize the energy. 
Spins in the (4,2) state are well aligned along the c axis. 
In the IC phase (represented by the (9,2) state) Sic. 
In the (5,2) state, moments are aligned mostly along the 
c axis but are slightly less oriented in (3,1) state until 
the field is close to H = 4. The spin reorientations at 
the boundaries (4,2)-IC and IC-(5,2) thus each represent 
spin-flop transitions. All transitions are first order. In 
the IC phases, planar components s x _L s y have equal 
magnitude and the biquadratic antisymmetric exchange 
coupling stabilizes a helically polarized spin structure. It 
is only in the IC state that an induced electric polariza- 
tion P _L c occurs^ 

Figures 5-7 illustrate the spin structures in the trian- 
gular layers at representative field values for each of the 
ordered states. The high-field (3,1) phase (HI) has a 
similar structure for each trilayer. In the zero field (4,2) 
state (TTi!) spins on subsequent trilayers alternate in di- 
rection. In the (5,2) phase, each of the six triangular 
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FIG. 3: Lower energy phases (L,M) as a function of magnetic 
field from the model Hamiltonian (3) using parameters Ji = 
1,J 2 = 0.33, J s = 0.3, J' = 0.4, Ci = 0.3, C = 0.1, = 
1.0, Gi = 0.06, G' = 0.02 with 3% axial anisotropy. 





FIG. 5: Schematic representation of three triangular layers 
(distinguished by color and font size) of the (3,1) and (4,2) 
phases where + and — symbols denote spin vectors along the 
c-axis. In the (4,2) phase, all spin vectors are reversed in the 
subsequent trilayer (m — 2). 





FIG. 4: Site averaged values of (sf) 2 corresponding to the 
lowest energy states from Fig. 3. 

layers in the magnetic unit cell are of the form (TTTII) 
but the relation between spins on alternating trilayers 
m = 1 and m = 2 is more complicated, as found by 
Mekata et alM- In the IC phase represented by the (9,2) 
structure of Fig. 7, spins lie in the triangular planes and 
rotate from site to site with a non-uniform angle. 

Model results on the evolution of the magnetization 
M z = (l/N)J2i S* with applied field is shown in Fig. 
8 and may be compared with corresponding experimen- 
tal data from Ref. p|. Key features of the data are 
reproduced by the model such as M z = in the (4,2) 



FIG. 6: Schematic representation (as in Fig. 5) of spins on 
the two sets of stacked trilayers characterizing the phase (5,2) 
at H = 2.5. 



state, increasing magnetization with field strength in the 
IC phase, plateaus at 1/5 and 1/3 in the (5,2) and (3,1) 
states, respectively. Flatter plateau regions, as seen ex- 
perimentally, can result from the present model by in- 
creasing the strength of the anisotropy or biquadratic 
symmetric exchange (see Discussion below). Note also 
that M z begins to deviate from the 1/3 plateau at fields 
above H C 4 ~ 4.5. This represents the transition to the 
canted (3,1) state proposed in Ref.[l|. The five critical 
fields resulting from the present analysis may be com- 
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FIG. 7: An illustration of spin vectors (in the plane) on part of 
the lattice representing (9,2) magnetic structure at H — 1.25. 
All spins are reversed in the subsequent tri-layer (m = 2). 

pared with the experimental values 7 T, 13 T, 20 T, 34 
T, and 70 T by multiplying the model field H by a fac- 
tor (J\S / gnu) — 5.6T. The resulting model critical-field 
values, 4.2 T, 9.8 T, 17 T, 28 T and 75 T, are in fair 
agreement with those deduced from the data, especially 
since no particular effort was made to adjust parameters 
for this purpose. 

V. DISCUSSION AND CONCLUSIONS 

A crude estimate of the contribution to biquadratic 
antisymmetric exchange interaction due to magnetoelec- 
tric coupling can be found using data on the dielectric 
constant e ~ 20 and polarization P ~ 400/iC/m 2 from 
Ref. 1]. With three Fe 3+ ions per unit cell&, the 
corresponding energy per magnetic ion is a very small 
Ep/N ~ P 2 /(e e) ~ 0.003K. Equating this to the an- 
tisymmetric exchange energy Eqp ~ C 2 /(8Ap), yields 
the crude estimate C\ ~ 0.1, somewhat smaller than the 
value used in the present analysis C\ — 0.3. Recall, how- 
ever, that the contribution Hcp to the spin Hamilto- 
nian can be deduced from general symmetry considera- 
tions and can in principle have a variety of microscopic 
origins, possibly stronger than from the magnetoelectric 
effect, such as higher order hopping analogous to sym- 
metric biquadratic exchange . 22 ! 29 (Note that an analysis 
of the present model with both stronger antisymmetric 
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FIG. 8: Site averaged values of (sf) 2 corresponding to the 
lowest energy states from Fig. 3. 



exchange and stronger anisotropy leads to an increase 
in all the transition fields as well flatter magnetization 
plateaus). Spin- flop transitions occur when the Zeeman 
energy Ez = — S • H equals the anisotropy energy Ea 
of a magnetic system, AE = E z — Ea = 0. Even very 
small perturbations can affect the system in the tran- 
sition region, particularly for highly frustrated systems. 
In the present case, the small antisymmetric biquadratic 
exchange interaction stabilizes the IC phase. This point 
is also demonstrated by considering an analysis of the 
model Hamiltonian in the absence of the term TLcp- The 
evolution of lowest energy states with C\ = C — is 
shown in Fig. 9, with all other parameters set as above. 
The IC phase as characterized above is never stabilized, 
even when considering a wider range of values for the 
other model parameters. 

Preliminary analysis of a Landau-type free energjsi^ 
based on a molecular-field treatment of the Hamiltonian 
(3) indicates that the proposed model captures essential 
features of the magnetic-field - temperature phase dia- 
gram. In this formalism, the free energy is written as 
a functional of the spin density s(r) = m + Se lC *' r + 
g* e — iQ r w j iere m j s th. e uniform magnetization and S 
and Q represent polarization and wave vectors, respec- 
tively, of the modulated spin structure. In zero applied 
field, the observed sequence of transitions from paramag- 
netic to linearly polarized IC to P4 states is reproduced. 
The stability of the linearly polarized P4 phase is as- 
sociated with fourth-order zero-field Umklapp terms of 
the form [S • S] 2 A4q i g, where G is a reciprocal lattice 
vector. Similarly, sixth-order terms [m • S][S ■ S] 2 A 5 q : g 
enhance the stability of the field-induced linearly polar- 
ized P5 state. The higher field period-3 phase occurs due 
to terms of the form [m • S][S ■ S]A3Q g- The helically 
polarized IC phase is characterized by S _L S*— is a re- 
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FIG. 9: Lower energy phases as in Fig. 3 calculated with no 
biquadratic antisymmetric exchange. 

suit of contributions arising from TLcp- Detailed results 



will be reported elsewhere. 

The present work demonstrates that the complex se- 
quence of ordered phases observed in CuFe02 stabilized 
by increasing magnetic field strength is a consequence of 
a high degree of frustration due to a multitude of com- 
peting interactions. The usual basal-plane period-3 spin 
structures associated with the STAF are destabilized by 
longer-ranged intra-plane as well as inter-plane quadratic 
exchange interactions. Weak axial anisotropy is found to 
be enhanced by biquadratic symmetric exchange, likely 
due to magnetoelastic coupling. Magnetoelectric interac- 
tions through spin-orbit coupling provides for one mech- 
anism that gives rise to a new type of biquadratic an- 
tisymmetric exchange term which serves to stabilize the 
helically polarized spin-flop phase. The model Hamilto- 
nian proposed here will serve as the foundation for future 
work focused on magnetoelastic effects, spin excitations 15 
and finite-temperature effects. Many of the features of 
the present model are likely relevant to other recently 
examined magnetoelectric STAF's^ 
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the Natural Science and Engineering Research Council of 
Canada (NSERC) and the Atlantic Computational Ex- 
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